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We provide a consistent description of the kinetic equation with triangle anomaly which is com- 
patible with the entropy principle of the second law of thermodynamics and the charge/energy- 
momentum conservation equations. In general an anomalous source term is necessary to ensure 
that the equations for the charge and energy-momentum conservation are satisfied and that the 
correction terms of distribution functions are compatible to these equations. The constraining equa- 
tions from the entropy principle are derived for the anomaly-induced leading order corrections to 
the particle distribution functions. The correction terms can be determined for minimum number 
of unknown coefficients in one charge and two charge cases by solving the constraining equations. 



The experiments at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) 
provide several pieces of evidence that the quark gluon plasma (QGP) may have been generated and is strongly 
coupled or nearly a perfect fluid (called sQGP), contrary to the conventional picture of the QGP as a weakly 
interacting gas of quarks and gluons (for reviews, see e.g. [H-[|[). Relativistic hydrodynamics is a useful tool 
to describe the space-time evolution of the fireball formed in heavy ion collisions. The RHIC data for collective 
flows have been well described by the ideal and dissipative hydrodynamics [8l-ll5l|. 

The correspondence of relativistic hydrodynamics to charged black-branes was investigated by AdS /CFT duality 
[Sill- A new term associated with the axial anomalies was found in the first order dissipative hydrodynamics 
(see, e.g., Ref. [IH about holographic hydrodynamics with multiple/non-Abelian symmetries, or Ref. pj| in Sakai- 
Sugimoto model). Recently the new term has been derived in hydrodynamics with a triangle anomaly [20]. A 
similar result was also obtained in microscopic theory of the superfluid |21l. This problem is closely related to 
the so-called Chiral Magnetic Effect (CME) in heavy ion collisions [22l - |25| |. When two energetic nuclei pass each 
other a strong magnetic field up to 10 18 G is formed, which breaks local parity via axial anomaly. This effect may 
be observed through charge separation. Hydrodynamics in an external background field can be used to pin down 
the CME in real time simulation. However the anomalous term in the charge current breaks the second law of 
thermodynamics unless new terms of vorticity and magnetic field are introduced in the charge and entropy currents 



In this paper, we try to provide a consistent description of the kinetic equation with a triangle anomaly. We 
will derive the kinetic equation to the next to leading order as well as the leading order correction to the particle 
distribution function arising from anomaly. These results are compatible with the entropy principle of the second 
law of thermodynamics and the charge/energy-momentum conservation equations. 

This paper is organized as follows. In Sec. UH we will derive constraining equations from the entropy principle 
for the correction terms in distribution functions in the most simple case with one charge and one particle species 
(without anti-particles). In Sec. IIII1 we will show that an anomaly source term is necessary in general to ensure 
that the equations for the charge and energy-momentum conservation are satisfied and that the correction terms of 
distribution functions are compatible to these equations. In Sec. HVI andfVl we will solve the constraining equations 
to obtain the correction terms of distribution functions in the one-charge (with anti-particles) and two-charges 
cases. Finally we summarize and make conclusions in Sec. I VII We adopt the convention for the metric tensor 
<^ = diag(+, -,--). 

II. CONSTRAINING DISTRIBUTION FUNCTION WITH ANOMALY COMPATIBLE TO SECOND 

LAW OF THERMODYNAMICS 

In this and the next sections we will consider the most simple case with one charge and one particle species 
(without anti-particles). The relativistic Boltzmann equation for the on-shell phase-space distribution f(x,p) in a 
background electromagnetic field ~ d il A v — d v A^ is given by, 
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where the charge of the particle is Q = ±1. Here p denotes the on-shell 4- momentum satisfying p 2 — m 2 where 
m is the particle mass. We note that C[f] contains a normal collision term Co[f] and a source term from anomaly 
Ca[J], C[f] = Co[f] + CaI/]- We assume that Ca[J] is at most of the first order, a small quantity. The necessity for 
the source term is to make the charge conservation equation hold, 

d^f = -CE^B^ = -CE ■ B. (2) 

Here j M is the charge current and E^ = u v F IJ ' v and B^ = he^a^u" F a ^ are electric and magnetic field vectors 
respectively, where it M is the fluid velocity and e^ va p = —e^ va ^ = —1,1 for the order of Lorentz indices {iivafi) 
is an even/odd permutation of (0123). However, the presence of the source term should not influence the energy 
momentum conservation, 

d^ v = F v »fr. (3) 
One can verify that the equilibrium solution of the distribution function, 



exp[/3u M (p^ — QF^Xy) — j3Qpo\ — e 

satisfies the collisionless Boltzmann equation (fTJ) in an external field for constant (3 — 1/T (T is the local tempera- 
ture), u M and fig (local chemical potential without electromagnetic field). Here e = 0, ±1 for Boltzmann, Bose and 
Fermi distributions respectively. When /3, and /xo are not constants but functions of space-time, the Boltzmann 
equation (fTJ) is not satisfied automatically. Note that we can absorb —Qx^u^F^ = Qx ■ E into fiQ so that /o has 
the form of an equilibrium distribution function, 



M*>P) = e ( u .p-J,)/r_ e ' ( 5 ) 

where /x = fiQ — x ■ E. 

We assume that the distribution function / in presence of an anomaly is a solution of the Boltzmann equation 
with collision terms in Eq. ([1]), where j3, it M and fi are functions of space-time. Generally f(x,p) can be written in 
the following form, 

f( X >P) = e (u- P -Q»)/l +x ( x ,p) _ ~ = /o( X 'P) + flix,p), (6) 

where fo(x,p) is given in Eq. (|5|) and f\{x,p) is the first order deviation from it, 



fl(x,p) = -f (x,p)[l + ef (x,p)}x(x,p). (7) 

It is known that a magnetic field is closely related to a charge rotation characterized by vorticity. So we introduce 
into the distribution function terms associated with the vorticity-induced current uj^ — ^^ va fiU l 'd a xiP and the 
magnetic field 4- vector B^ which are assumed to be of the first order, which provide a leading order correction to 
the particle distribution function. For simplicity we will neglect viscous and diffusive effects throughout the paper, 
then the ordinary form in the current scheme for x{ x iP) reads, 



X{x,p) = \{p)p-uj + \ B {p)p-B, (8) 

where \{p) and Xb(p) are functions of /i, T, and u ■ p and have mass dimension —2 and —3 respectively. We will 
show that X(p) and \b(p) must depend on momentum otherwise they will contradict the entropy principle from 
the second law of thermodynamics. 

Using Eq. ^ we can decompose the charge and entropy currents and the stress tensor into equilibrium values 
and the leading order (first order) corrections as ^ = + jf , 5" = + and T» v = Tg v + with 

#i(aO = Q f[dpWfo,i{x,p), 

S$(x) = - J[dp}p^j(f ), 

S?(x) = - J [dp}p^'(f )h, 

Tft(x) = j[dp]p»p v f 0A {x,p), (9) 
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where we have defined [dp] = d g ( 27r ^ ( p M . p) (d g is the degeneracy factor), tp(f ) = / ln(/ ) - e(l + e/ )ln(l + e/ ) 
and vp'(fo) = ln[/o/(l + e/o)] = — (u ■ p — Q/j,)/T. Inserting /o into the above formula, we obtain the charge and 
entropy currents and the stress tensor in equilibrium, j'q = nu^, Sq — su^ and Ttf v = (e + P)u^u 1 ' — Pg^ u , with 
the energy density s, the pressure P, the particle number density n and the entropy density s = (e + P — n^i)/T. 
Using Eqs. (|7l8l9p , we obtain 



T? v = DT(u^ui v +u v u)^) + D B T(u^B v + u"B»), 



■^u" + Z B B») + (£>w" + DbB"), 



(10) 



where 



D 
D B 



-QJ£i = ^Qj [dp][(p ■ u) z - m J ]f (l + e/ )A(p), 
-Q4i B =\Qj [dp][(P ■ uf - m 2 ]/ (l + ef )\ B (p), 



_31 = 

T ~ 3T 



[dp][(p ■ uf - m 2 ](p • u)/ (l + e/ )A(p), 



J 31 



T 3T 



— / [dp][(p • u) 2 - m 2 ](p ■ «)/ (l + e/ )A B (p). 



(11) 



On the other hand, £, £b, -D and Db as functions of \x and T can be determined by the second law of thermo- 
dynamics or the entropy principle together with Eq. <[2][3j) . We find that d^(su^ + S^) cannot be positive definite 
unless we make a shift to introduce a new entropy current S 11 as follows, 



= su" + S[ l - (Dui» + D B B») = su^ - ^ + i B B» 
= ^(Pu»-^ + u x T x »)-(D^ + D B B>*) 



(12) 



We will use the thermodynamic relation 



d ll {Pu> 1 )=tid l JI-T^d ll u x 



(13) 



and the identities 



u^u^duUJx 



M e + P 



(nw^E^+cj^P), 



d^B^ = 2u p E p - —— 



(nB x E x + B»d p P), 



(14) 



to evaluate c^S^. We have used the shorthand notation /i = fi/T in Eq. JT3|). Following the same procedure as in 
Ref. [20], we obtain 



d p S p = uj» 



2-D „ 

7Tp- d » p 



+B" 
+E ■ oj 
+E-B 



Cl s a M 7Z - d^D B + -^d„.P 



D B 

7Tp°" ! 



1 ™ 2nD 

ss ' -2D f 



e + P 



1 tSS , r PA , nD B 



(15) 
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where we have defined 



fctF-T^-fc- (">) 

For the constraint d^S 11 > to hold, we impose that all quantities inside the square brackets should vanish. We 
finally obtain 

1 u 3 1 u 2 sTu 2 sTu , , 

Using Eqs. (|16|17| . one can verify that the values of £ ss and £g S are identical to Ref. (20j . The difference between 
our values in Eq. IT7|) and those in Ref. [20| arises from the fact that we do not use the Landau frame while the 
authors of Ref. [20| do. By equating Eq. (fTTj) and (fl7j) . we obtain equations for A and Xb, 

QJ& = J A = — DT, QJ 2 A f = -£ B , J 3 A f = -i? fl r. (18) 

Equation (I18p forms a complete set of constraints for A and Xb- We note that A and Xb must depend on momentum 
in general. If A and Xb are constants, we would obtain 

X = __*_ = Q___, (19) 
DT D B T *J 31 ' 1 7 

which contradict Eq. (|17[) from the entropy principle. 
We can expand X(p) and Xb{p) in powers of u ■ p, 

A(p)=£> («•!>)*> A B (p) = (20) 

i=0 i=0 

So we obtain the following expressions 

Jnl = ^^hJi+n.l, J„l = Jj+n,U (21) 

i=0 i=0 

for n = 2, 3. Here the functions J ng are integrals defined in Ref. [7l [TTj. 

Jn q = (-iyj^ I ^^[(-^--^(-^/oCl + e/o), (22) 

Using Eqs. (|2QI2ip in Eq. (fT5)) . we can constrain the coefficients A^ and Xf . If we expand both X(p) and Xb(p) 
to the first power of u ■ p, we can completely fix the coefficients A0.1 and A^ from Eq. (fT5|) since we have two 
equations for Aq,i and two for Xq X , 



whose solutions to Aq 1 are 



( QJ21 QJ31 \(*o\ ( ~t 

\ j 3 i Jax ) v A i y V - DT 



X \ _ 1 f -CJu + DTQJ31 



Ai / Q(J 2 i J 4 i - V ^ J 3i - DT QJ2i 



(23) 



(24) 



The equations and solutions for A^ are in the same form as Eqs. (|23|24j) with replacements Ao,i — > A^, £ — > £b 
and D — > Db- 

For massless fermions and small /Z, the results are 

A « -CGi^j, Ai « CG 2 ^5, Aq 3 w — , Af w — , (25) 



where Gi and G 2 are two constants, Gi = 6075007r 2 C(5)/(d ff G ) and G 2 = 12607r 6 /(d g G ) with G = 
455625C(3)C(5) — 497r 8 . We notice that the D terms in Eq. (12"4"|) are negligible, so the solutions are proportional to 



III. COLLISION AND ANOMALOUS SOURCE TERMS 



In this section we will show that a general form of A(p) and Xb(p) are compatible to the charge and energy- 
momentum conservation equations ([2]) and ((3j) . We will also derive equations for the collision and anomalous source 
terms. For simplicity we consider the single charge case without anti-particles. 

The charge conservation equation ^ can be derived from the Boltzmann equation ([T]) as 



where have used the identity 



(2tt) 3 E p 

d 3 p 
(2ir) 3 E p 



d "P p » F K 
(2n) 3 E p P ^ v dp v 



{2tt) 3 E p 



C[f], 



(26) 



= - J d 4 p^-[9(po)S(p 2 -m 2 )p^F^f = Q. 



Then the momentum integral of the collision term must obey 

d 3 p 



(2w) 3 E p 



C[f] = -CE-B, 



(27) 



(28) 



so that Eq. @ can hold. The energy and momentum conservation equation Q can be derived from the Boltzmann 
equation (J} as 



d p T^ 



d 3 p 



(2tt) 3 E 1 



d 3 p 



wkr/ paFaP 



df 



dpp J {2^E p 



d 3 P v r \t\ 



(29) 



where we have used 



2E n 



p»p a F a p-^ = d 4 p9{p )S(p 2 -m 2 )p»p a F a0 



df_ 

dpp 







d 4 p^[6(po)6(p 2 - m 2 )p»p a F af3 ]f = F» a j 



Then we require that the collision term must satisfy 



d 3 p 



(2nfE p 



p v C[f\ = 0, 



so that Eq. © can hold. 

We can expand C[f] to the second order as 

C[f] = C [/o + fi + M + C A [f] « C x + C 2 
where we have used the property Co[/o] = 0, and defined 



(30) 



(31) 



(32) 



Ci = 

C 2 = 



dC 
df 

dC 
df 



f=fo 



H+Cax, 

. 1 d 2 C 
12 2 df 2 



f?+C A 2- 



(33) 



f=fo 



Note that the general form for the normal part of C\ is ^# f\ = H\(u-p)p-uj+H\ B [u-p)p-B. When inserting the 



f=fo 



distribution function ([5]) into the Boltzmann equation ([T]) and using Eq. (|32l) . we obtain the Boltzmann equations 



G 



to the first and second order, 



\ dx* 1 M dp v 



fa — Ci, 

h = c 2 . 



From Eqs. ()33|34|1 we can determine the anomalous source term of the first order, 

C A i = ~H\{u -p)p-u3- H\ B (u ■ p)p ■ B — / (1 + ef )p^d^[(u • p - fJ,o)/T]. 
By evaluating the left hand sides of Eq. (|35[) , we can fix C 2 as follows 

C 2 = /o(l + e/o){-[(l + 2e/o)Aa^ , (/o) + (5 Al A)]pV^ 
-[(1 + 2e/ )A B 9^'(/o) + {d ll X B )WB v 
dX 



(34) 
(35) 

(36) 



-/3(l + 2e/ )A + 
-13(1 + 2e/ )A B - 



d(u ■ p) 
d\ B 



d(u ■ p) 



p^E^Wv + Ap^F^w" 
p»p u E„B v + Xbt^F^B" 



(37) 



Taking momentum integrals for Eqs. (|34I35)1 . we obtain the divergences of charge currents to the first and second 
order, 



d^(x) = J [dp]Ci = J [dp]C A1 , 
M» = J [d P ]C 2 . 



where we have used the property, 



f [M h = I W] [Hx(u ■p)p-iu + H Xb (u ■ p)p ■ B] = 0. 

J aj f=h J 



(38) 
(39) 

(40) 



With Eq. (JQJ) and identities in Eq. (|T4j) . the left hand side of Eq. (|39j) can be evaluated as 



2n£ 
e + P 
2£ 



B ■ E 



e + P 



e + P 



d.P +5" d^B 



e + P 



zd u P 



(41) 



In order for Eq. @ to be satisfied, it is required that d^(x) must have the form 
d^{x) = -CB-E-d^(x) 



lj ■ E 2£s - 



2< 



e + P 

2e 



BE C 



e + P 



e + P 



8.P -F cUb 



e + P 



,d u P 



(42) 



We see that <9 M jji (x) is not vanishing but of the second order though it is superficially a first order quantity. This is 
very important otherwise it would lead to d^j^(x) = —CE ■ B and give rise to incompatible results for £, £b, D, D b 
with Ref. [23|. Therefore we should keep in mind that it is the full charge current that satisfies charge conservation 
equation ©. We note that the expression of dajn(x) in Eq. (|4"2")l also gives the momentum integral of the anomalous 
source term of the first order, 



[dp]C 



Al 



-u ■ E I 2& 



2< 



e + P 
2£ 



B-E C 



P 



d^P ) - B" ( d^ B 



e + P 
t 



d u P 
e + P 



(43) 
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For the energy and momentum conservation we obtain, 

d^-F^jl = J [dpWd, (44) 
fyZT - F^jl = J [dp]p»C 2 , (45) 

With Eqs. (|9ll4p . we can evaluate the left hand side of Eq. (|4"5")l as, 

g^Tf _ F »»ji = U j v [DTd-u + u-d(DT)] + B v [DBTd-u + u-d(D B T)} 

+ u«{. • E(2D B T -^-E. + w • [d{DT) - |f^P] 

+p.[9(P B T)--^9P]} 

+DT(u -du v +uj- du v ) + D B T(u ■ dB v + B ■ du u ) 

-F^^ + ZbB^ (46) 



where one can verify that each term in the right hand side is of second order. The left hand side of Eq. (|44j) is then 
given by 

W - F^jI = -(W - F^jl), (47) 

which is also a second order quantity though p v C\ = p v p^ {^j: — F^ a g^-^ /o is superficially of first order. One 

might question the validity of Eq. (fT4")l which follows d^T^" — F v ^j^ = 0, but it is not a problem here since this 
equation really holds at the first order or the leading order but not true at the second order. It is essential that 
the energy momentum equation Q hold for the full quantities and j M , and not for Tq\ and j®' 1 separately, 
otherwise the results would be contradictory to those of Ref. [20] following the entropy principle of the second law 
of thermodynamics. 

We now obtain the momentum integral of p v C\$ from Eqs. (|44I45|) with d^Ttf^ — F^^j®' 1 given by Eqs. (|46|47j) . 



IV. ONE CHARGE WITH PARTICLE/ANTI-PARTICLE 

We now consider one charge case but add anti-particles to the system. We will calculate X(p) and \b(j>)- Since 
there arc particles and anti-particles, we recover the index Q = ±1 in particle distribution function, / — > and 
/o,i — > foi- In Eq. (J9j) , summations over Q should be added. In Eq. (JTlJ we also have to add the index Q to 
J^x and Jni- J.„i — > an d Jni — ► J^i ^ ( n = 2,3), and add summations over Q into the formula of £, £g, 

D, D B - Inserting /q into Eq. we obtain the charge and entropy currents and the stress tensor in equilibrium, 
jg = fin'", Sq — su^ and Tq V = (e + P)u^u v — Pg^ u , with the particle number density n = J^q Q n Qi the energy 
density e = £q, the pressure P = Pq, and the entropy density s = sq = (e + P — n/x)/T. The solutions 
to £, £b, D, Db are the same as in Eq. (|17[) . Then Eq. (|T5)) is modified to 

tA,+ tA,- _ p tA,+ , tA,+ _ _Jjrp 

J 2l J 21 ~ Si j 31 ' .31 — ' 

J 21 B ' + - J 21 B ' = -£b, ^3i B,+ + ^3i B: =-D B T. (48) 

Eq. (j2"Tj) now becomes 

tA,Q _ \ tQ tA b ,q _ \p <B tQ /.qx 

i=0 i=0 

In the case of minimal number of coefficients we can completely fix the coefficients Ao,i by solving following system 
of equations, 



6J21 SJ 31 \f A \ _ / -C 
<jJ 3 i oJax \ Ai J I -DT 



(50) 



where we have used the shorthand notation, (5J n i = — J nl and oJ n \ = J^i + J nl . The solutions to Ao,i are 

(51) 



M — ( _ £( ctJ 4i) + DT(SJ 3 i) 
Ai J A V £( CTj 3i) - DT(5J 21 ) 
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where A = (<5J2i)(crJ4i) — (&J3i){8J3i). The equations and solutions for Xq X are in the same form as Eqs. ([50I51|) 
with replacements Ao,i — > A^ 1; £ — > £s and D — > Dg. 
For massless fermions and small Jt, we have, 

7 4 

o.hi « 9C(3)G, <rJ 31 * -^—GT, aJ 41 « 225C(5)GT 2 , 

15 

7tt 4 

(5J 2 i » tt 2 71G, (5J 3 i « 36C(3)7ZGT, <5J 41 w — tjGT 2 , 

O 

7T 2 

A « —G 2 G T 2 -p. (52) 
5 

where G = and G = -84tt 2 C(3) + 1125C(5). The solutions have very simple form, 

A ~ CGx-j^, Ai w -GG 2 ^j, Ajf = -A Af = (53) 

where we have used two constants, Gi = 675 ° 7r G ^ 5 ^ , G 2 = ■ We notice that the D terms in Eq. (|5"Tj) are 
negligible, so the solutions are proportional to £. Note that the quantity -^p ~ £t 2 is also small and we have 
dropped it, since n[i ps ^-T 4 71 2 and e + P = |e « d g ^-T 4 . 

V. WITH TWO CHARGES AND PARTICLE/ ANTIPARTICLE 

As an example for the case of two charges, we consider adding to the system the chirality or an axial U(l) 
charge to particles. Then there are two currents, one for each chirality, or equivalently, for the U(1)/Ua(^) charge. 
For simplicity we assume that there is an anomaly for the axial charge current but no anomaly for the charge 
one. There are distribution functions for right-hand and left-hand particles, f®(x,p) (a = R,L), with chemical 
potentials Hr.l = n ± (J>a- As an extension to Eq. ©, the corrections Xa( x ip) m fa( x >p) are now Xa(x,p) = 
XaP^uj^ + XaBP^B^. Instead of right-hand and left-hand quantities X a , we can equivalently use X = X R + Xl 
and Xa — Xr — Xl, where X — A, As, £, n, s, e, P, j M . The distribution functions fQ(x,p) satisfy two separate 
Boltzmann equations of the following form, 

p " {-L qf ^w) f « {x > p) = CaQ[f »' ] - (54) 

The U(l) /Ua{X) charge and entropy currents and the stress tensor in equilibrium are give by: j§ — = uau^, 

Sq = su^ and T$ v = (e + P)u^u v — Pg^ h ' , with the particle number density n = J2 a Q Q n aQ {a — R,L; Q = ±1), 
the energy density e = J2 a Q e aQ, the pressure P = J2 a Q p a.Q, and the entropy density s = J2 a Q s aQ = ( e + P ~ 

T,a n ^a)/T. 

Similar to Eqs. (|26l28p . the J7(1)/J7a(1) charge conservation equations can be derived as 

d»f(x) = J[dp]{C R , + -C R , + +C L , + -C L ,+)=0 1 (55) 
dMx) = J [dp]{C R .+ - C fl ,_ - C L .+ + C Lt+ ) = -CE»B^ (56) 
where we have used Eq. (|27[) . The energy and momentum conservation equation reads, 

d ^ v _ = J [dp]p ^c R + + Cr _ + C L .+ + C L ,_) = 0, (57) 

where we have used Eq. (|3U|) . 

Similar to Eq. pip , we can express £ a , £ a s, D and Dg in terms of A a and X a s as, 

t _ _( 7-A a ,+ t x °'-) _ dt — J Xr ' + _|_ j x n - 1 j x l,+ , tAl,- 
C|a — W21 ^21 )i 1J± — J Z1 ' J 31 ' J 31 ' J 31 > 

e aB = -(4° b - + - 4r n, -d b t = 4r ,+ + 4* B '~ + 4r ,+ + & B, ~, 
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for a = R,L. We have the following first order corrections, 



Jl 



Jai 



a=R,L 



im ~ Hx = (& ~ + (Zrb - ZlbW = + UbB», 

DT (/w" + u'w") + D B T (u^B v + u y B") , 



E m (£a 



a=R,L 



T 



= - E 



i=n\i\\,A 



T 



" +& B B») + {Dw»+D B B> 1 ). 



(59) 



It can be verified that the entropy current in Eq. (|59|) cannot satisfy d^S^ > unless C = 0. In order to ensure 
the positivity of d^S^- in presence of an anomaly, one would have to substract the vector = Duj^ + D B B^ from 
S*. With U{1) and U A {1) charges, we have 



S" 



S" -Q" 



E TlaiCa^ + ZaBB*) 



a=R,L 



(60) 



a=R,L 



In the same way as in Sec. HU the divergence of the entropy current can be evaluated as, 



E d ^a 



a=R.L 



n{TD 
e + P 



2D 



-B" 



E d »Va 
a=R,L 



_ ( n{TD B 



+ E-LJ 



+E ■ B 



E 

a=R,L 

E 

a=R,L 



e + P 

n a D _ 
i + P T 



n a D B _ £aB_ 
e + P T 



-Lb) - d^D B 



D, 



2nD 
e + P 



-2D T 



-OaP 



u a nD B 

c ^F + 7Tp 



(61) 



following Eqs. (|55I(|57p). By imposing all quantities inside the square brackets to vanish we can solve £ a , £ a B, D, 
D B as follows, 



£ B = CfJLA 1 - 



2e + P 
2nfi 



Cab = Cull - 



2n A /j, A \ 
e + P J 



e + P , 

DT = -Cha^ 2 , D B T = -ChaH 
It is interesting to observe that for small /i, fiA, we obtain very simple and symmetric solutions, 

£ k 2Chha, £,a ~ Cy 2 , £ B w C/j, A , Cab ~ C/j,, 
DT = -Cha^ 2 , D B T = -CfiAfJ>- 



(62) 



(63) 



which is identical to the result of Ref. |23l . |25j. Here we have assumed that all integral constants are vanishing. 
Note that Eq. (|63p is the result of the entropy principle in the hydrodynamic approach which was also obtained in 
Ref. gHH. 

Equivalently we can use Cr = (£+£a)/2 and Cl = (£~£a)/2 to determine A a , X aB (a = _R, L) via solving a system 
of equations (|58|) . where the first/second line (each has three equations) is for \r,l/^rb,lb- For minimal number 
of coefficients we can determine the values of these coefficients completely, Xr l an d X bb> l B can be expanded to 
the zeroth or first power of u ■ p. For example, if we expand Xr to the zeroth power, then we have to expand Xl 
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to the first power, and vice versa. Suppose we take the former case, A^ = A^ and \l = ^lo + Xli{u ■ p), we can 
solve Xr,l as 

A_ro = —£r 



sj 2 r 



Alo = ~Zl{<tJ&) + 1 (DT Zn^jP) (SJ 3 \), 

Ali - ^L(aJt) -\{dT- & ^|) (^2 L i), (64) 

where A = (<5^2i)( cr ^ii) — (<^3i )( CT ^3i )• The solutions to Xrb,lb take the same form as above with replacements 
A a — > A a s, £ a — > and D — > if we assume the same expansion as Xrx'- Xrb — Xrb,o and Xlb = 
Xlb,o + X LB .i(u -p). 

For massless fermions and small ~p R L (or equivalently small /x and Ha), we obtain 

<tJ& l « 9C(3)G, *jg L « ^GT, aj£ L « 225C(5)GT 2 , 

<5j£ L « Tr 2 ~p R L G, 5J£ L » 36C(3)7Z«, L GT, 5J 4 ^ L » L GT 2 , 

A « ^G 2 G T 2 m l , (65) 
where G and Go are the same as in the former section. Then the solutions to \r,l are 

3G (i (i + 2\i A 3C fi f u + 2fi A n — 2/j, A 
aro ~ — j~7Pi ; > a lo ~ ~r ^4 ; h Cr 5 



A 



LI 



d g T 3 /x + ha dg T 3 \ h + ha li- ha 

lAit^C /x HHa 



3G 1 3G 1 307t 2 G M Ma 

Xrb.o « -^^3, A LB , « — ^, Alb,i « ^7^' (66) 

where G 4 = — ^-847T 2 C(3) and G 5 = ^1125^(5). The coefficient ratios of X ai /X aB .i (x = 0, 1) are proportional to 
H times dimensionless factors, 

Aflo H + 2 MA A L0 / ~ H + 2/xa , _ M - 2 MA 
/x , w /x G 4 h G. 



AfiB,0 M + A* A X L b,q \ H + HA li- HA 

Ali 7 2 T 



2 



- W>2 2- ( 6? ) 

Als.i 15 h 2 - H A 

Note that Xlb,i "C TXlb.o ~ TXrb.o, so both Xrr{p) and Xbl{p) can be constants at small 77 and /Z4 limit. 
This property is quite different from the one-charge case in which Xb(p) must have momentum dependence in order 
to comply with the entropy principle. 



VI. DISCUSSIONS AND CONCLUSIONS 



We have shown that induced terms related to the vorticity and magnetic field in the charge and entropy currents 
from a triangle anomaly can be derived in kinetic theory by introducing correction terms to the phase space 
distribution function at the first order. We demonstrated that the anomalous source terms are necessary to ensure 
that the equations for the charge and energy-momentum conservation are satisfied and that the correction terms of 
distribution functions are compatible to these equations. 

As examples for the correction terms of distribution functions, we focus on the massless fermionic system in three 
cases for small /x/T, with one charge [L/(l)] and one particle species (without anti-particles), with one charge [J7(l)] 
and particles/anti-particles, and with two charges [1/(1) x [/^(l)]. In the latter two cases, the coefficients for u> and 
B terms in distribution functions are found to be proportional to G/i/T 3 and G/T 3 respectively. In the two-charges 
case the coefficients can be constants or independent of momentum, such a property is impossible for the one-charge 
case since it is not allowed by the entropy principle. 
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In the two-charges case, we assumed that there is an anomaly for the axial charge current but no anomaly for 
the charge one. The coefficients of correction terms for the charge/axial-charge currents and energy-momentum 
tensor have a very simple and symmetric form at small fi/T and ha/T limit: £ w 2C^,^a, ^a ~ C/i 2 , £b ~ C\ia, 
£,ab ~ Cfi, DT = —ChaH 2 , and DbT = —CfiAH- This means that similar to the CME an axial anomaly can 
induce a residual charge current which is proportional to the magnetic field and the axial chemical potential (23| . 

We have a few comments about our results. In our evaluation of the correction terms of distribution functions, 
we have assumed that X(p) and Xb{p) are identical to particles and anti-particles. Alternatively we can assume 
that they have an opposite sign for particles to anti-particles. We can not tell which case is correct due to lack of 
deeper knowledge about these anomalous term at a microscopic level. Similarly we have assumed that X{p) and 
Xb{p) have different values for right-handed particles from left-handed ones. One can also assume that they have 
the same or opposite values for right-handed and left-handed particles. In the current framework one can not tell 
which is correct. Such a situation is like what happens in an effective theory when many effective candidates point 
to a unique microscopic theory. We also note that the solutions to X(p) and Xb(p) given in this paper are for the 
cases where the number of unknown coefficients in X(p) and Xb(p) is equal to that of constraining equations. It is 
possible that X(p) and Xb{p) can be expanded to higher powers of (p ■ u) and then have larger number of unknown 
coefficients than that of constraining equations. In this case the constraining equations just provide constraints for 
X(p) and Xs{p) from the second law of thermodynamics. 
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